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CN ■ Abstract. Given independent normally distributed points A, B, C, D in 

Euclidean 3-space, let Q denote the plane determined by A, B, C and D denote 

the orthogonal projection of D onto Q. The probability that the tetrahedron 

ABCD is acute remains intractible. We make some small progress in resolving 

this issue. Let T denote the convex cone in Q containing all linear combinations 

A + r {B — A) + s (C — A) for nonnegative r, s. We compute the probability 

that D falls in (B + C) - T to be 0.681..., but the probability that D falls in T 

PLh ■ to be 0.683.... The intersection of these two cones is a parallelogram in Q twice 

r^ the area of the triangle ABC. Among other issues, we mention the distribution 

"t3 ' of random solid angles and sums of these. 

Let oi, 02, a 3 , &i, fo 2 , &3, ci, C2, C3, d\, di-, d^ be independent normally distributed 
random variables with mean and variance 1. The points 
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A = o 2 , B = 62 , C = 02 , £) 





constitute the vertices of a tetrahedron in Euclidean 3-space. The tetrahedron ABCD 
is acute if each of its six internal dihedral angles are less than tc/2. It is known [1] 
that ABCD is acute if and only if the orthogonal projection of each vertex onto the 
plane of the opposite face lies within that face. Neither characterization suggests 
an easy approach to finding the probability that random ABCD is acute. We will 
examine a variation of the latter characterization, focusing on the point D. Let Q 
denote the plane determined by A, B, C and D denote the orthogonal projection of D 
onto Q. While calculating the probability that D falls within the triangle ABC seems 
difficult, we succeed in computing the probability that D falls in either of two convex 
cones in Q containing ABC. In fact, the cones contain the parallelogram with vertices 
A, B, C, —A + B + C; the fourth vertex is clearly the vector sum (B — A) + (C — A) 
displaced so that it emanates from the point A. Moreover, the intersection of the cones 
gives precisely the parallelogram. Further research might uncover other probabilities 
for related cones in Q, and thus the inclusion-exclusion principle might yield the 
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probability that D falls within the parallelogram. This task is best left to someone 
else (!) but we hope that our work provides some inspiration along the way. 

In section 1, we review a well-known proof that planar triangles ABC are acute 
with probability 1/4. Returning to 3-space, we compute in section 2 the probability 
that D falls in the cone 

{-A + B + C) - {r (B - A) + s (C - A) : r > 0, s > 0} 

using the Krishnaiah bivariate F-ratio distribution [21 El II]- This is the cone with 
vertex —A + B + C and outgoing edges parallel to vectors A — B, A — C. In the 
following section, we compute the probability that D falls in the cone 

T = A + {r (B - A) + s (C - A) : r > 0, s > 0} 

using a forgotten result of Miller's [5] , extended via convolution. This is the cone with 
vertex A and outgoing edges parallel to vectors B — A, C — A. The two probabilities 
are numerically close but not equal (at least not in our view). A rigorous proof of 
such an inequality is open. 

1. Triangles and Tetrahedra 
For the moment, 



A = W1 , B = r 1 , C 






are vertices of a random Gaussian triangle in Euclidean 2-space. Let Q denote the 
line determined by A, B and C be analogous to before. Also, let a,/3, 7 denote the 
angles at A,B,C respectively. The orthogonal projection of w — C — A onto the 
subspace spanned by v = B — A is 

y.w (B-A)-(C-A) 

v-v V ~ (B-A)-(B-A) {B A} 

and clearly C falls between A and B on Q if and only if 

{B -A).(C-A) 

(B-A).(B-A) 

Let 

^ P((B -,4).(c-.)>oHP( g:^g:^ >o). 
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then 



T ( (B-A).(C-A) ; 
F \(B-A).(B-A) / 



P ((B -A)-{C-A)<{B-A)-{B- A)) 

P((B-A).[(C-A)-(B-A)}<0) 

P((B-A)-(C-B)<0) 

p((A-B).(C-B)>0)=p 



by symmetry, hence 



P (°< {1~-A)'- { (B-A) <1 ) =P - {1 - P)=2P - 1 - 

To compute p, note that (B — A) ■ (C — A) is equal to 



3 / 2 , 5 C7 



3' 
1 



-A 



B C 




B C N 

72 + 72, 



_ 2 _ _ -2 
r)X2 r)X2; 



a linear combination of independent chi-square distributed variables (each with 2 
degrees of freedom) . Therefore 



/' = P ( 2*2 - 2*2 > 



KM 



ro 2 > — 

' 3 



3 
4 



by use of the F-ratio distribution with (2, 2) degrees of freedom. It follows that 
C falls between A and B on Q with probability 2(3/4) — 1 = 1/2. Equivalently, 
P(a > 7r/2) = 1/4 because cos(a) < precisely when (B — A) ■ (C — A) < 0. Since 
at most one angle of a triangle can be obtuse, we deduce that 



P (triangle ABC is acute) 



1 — P (triangle ABC is obtuse) 

1 - (p(a > 7r/2) + P(/3 > vr/2) + P( 7 > vr/2)) 
1 - 3/4 = 1/4 



as was to be shown. Our proof imitates Eisenberg & Sullivan's [6] approach, although 
we avoid angles a,(3,j until the last step, using C instead. Portnoy's [7] argument 
employs triangle medians rather than orthogonal projections. 
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Let us now return to random Gaussian tetrahedra in 3-space. Our argument is 
similar but more complicated. The orthogonal projection of w = D — A onto the 
subspace spanned by u = B — A, v = C — A is 

8 ' * a + ttv = ( ^-^^- a I(b - a) + { c,- a ); {d - a )ac - A) 



u ■ u 



V ■ V 



(B-A)-(B-A) K ' (C-A)-(C-A) 

and clearly D falls within the desired parallelogram on Q if and only if 

(B-A)-(D-A) , (C-A)-(D- A) 
< \- ^4^ Jr < 1 and < \- ^4^ Jr < 1. 



(B-A)-(B-A) 

This is the same as requiring that 

(A-B)-(D-B)> 0, 
(A-C)-{D-C)> 0, 



(C-A)-(C-A) 



(B-A)-(D-A)> 0, 
(C-A)-(D-A) >0. 



Each product is of the form (3/2)%| — (1/2)%| and jointly they give rise to corre- 
sponding F 3j3 ratios: 



A 2 D 

76 V3 + 7q 




/2 ^B_ ^D_ 

3 ^6 + v 7 ^ 



5 D 

72 72 



> 



/2 , C 


D 

76 




C D 

~V2 72 


2 



> 



A computation of the probability that all four inequalities hold simultaneously does 
not seem possible. We note that the two expressions in the left-hand column possess 
the same denominator, which is essential for the next section. An entirely different 
approach will be needed for the two expressions in the right-hand column. Finally, 
the two left-hand inequalities are true if and only if the orthogonal projection of w is 
ru+ sv for some r < 1, s < 1; this translates into D falling in the cone 

A + {(1 - r) (B - A) + (1 - s) (C - A) : r > 0, s > 0} = (B + C) - Y. 

Likewise, the two right-hand inequalities are true if and only if the orthogonal pro- 
jection of w is ru + sv for some r > 0, s > 0; this translates immediately into D 
falling in the cone I\ 
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2. Krishnaiah Bivariate F- Ratio Distribution 
Let X be an n x 2 matrix of n independent random row 2- vectors, each distributed 
according to N(0, E), where 



-1 < p< 1. 



Let &ij = Xi ■ Xj, where X/. is the k th column of X for k = 1, 2. Also let f = Y ■ Y, 
where Y is a random column m-vector satisfying Y ~ N(0, 1) and independent of 
X. Then the joint distribution of (man)/(nf), (ma22)/{nf) is a bivariate F-ratio 
distribution with (n, m) degrees of freedom and with S as the associated covariance 
matrix (for X). These conditions are clearly met in our case, for which n = m = 3, 



/ __^L A. \ 

a 2 4 

«3 ds_ 

V >/2 V2 / 



X= ^ 





y 



'2 4 

3 C3 + 7!/ 



and 



It follows that j3] 



Cov 






Varte + ^U- 



/2 4 a fc /2 4 

3 h + 7E-T6'i~3 Ck + 7E 



1 



pf^^and^ 2 -^ 



nr 



nr 



(1 - p 2 )"/ 2 ~ p 2fc r(n + m/2 + 2ifc) 

2^ rrrr7 /» . ,x A fc 



where 



A, 



r(m/2)r(n/2) fc f A;!r(n/2 + A;) 

^n/2+fc-l 



77 77 



;i + X + yY+^k+m/2 



dxdy. 



and 7] = (n£) / (m(l — p 2 )). Each A/, can be evaluated symbolically and the series 
appears to converge fairly quickly. For our case, £ = 1/3, hence r\ = 3/8 and the 
desired probability is 0.6810669069.... 
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3. Miller Bivariate Density for Products 
We want the probability that both (B - A) ■ (D - A) > and (C - A) ■ (D - A) > 0. 
For now, we consider a simpler scenario in which a = A,b = B,c = C,d = D are 
scalars. 

Let X be a random p- vector and y be a random scalar. Assume that (X, y) ~ 
N(0, E) and that the {jp + 1) X {jp + 1) covariance matrix £ has inverse 

where Q is p x p, v is p x 1 and u is a scalar. Let Z = y X. Then the joint density 
of Z is 



ydetg-i) / w n(p-d/4 rr^o^ 

(27r) 0H-i)/ 2 (^^J ex P(-^ ^)^(p-D/ 2 (v« ^J 

where K( p _i)/ 2 (6 I ) is the modified Bessel function of the second kind. In our case, 
p = 2, 

(( xi\\ 1 1 b-a\\ / 2 1 1 

£ = Cov x 2 = Cov c - a = 1 2 1 

hence 

11 1 

^dettE- 1 ) = -, -u'Z = - (zi + z 2 ) , Z'QZ = - {3zf - 2z x z 2 + 3z 2 2 ) . 
Also __ 

Kl/2{e) = V 2 ^T" 

thus the density simplifies to 

1 exp ( - ( zj + z 2 - vW 32 ! - 2zi 2 2 + 3zf ) ) 



271 



oZi — 2Z\Z 2 ~4~ OZ 2 



In fact, we wish to compute the probability that the sum of three independent 
copies of Z has both components > 0. One way to do this is to evaluate the sextuple 
integral: 

00 00 00 00 00 00 

////// f( Zl ~ Ul ~ Vl ' z 2-U2-v 2 )f(ui,u 2 )f(vi, v 2 ) duidu 2 dvidv 2 dzidz 2 = 0.683. 

—00 —00 —00 —00 
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obtained via repeated convolution of / with itself. Another way uses the Fourier 
transform 



F(wx,w 2 ) 



1 



\J(w\ - i)(3iui + i) + (w 2 - i)(3w 2 + i) + 2witu 2 - 1 



which, when cubed, maps back to a remarkably simple density function. Higher 
precision is now possible: 



o o 



4-v/3' 



exp 



7T 



z\ + z 2 — \/2>J?>zl 



2z±z 2 -\- oZ 2 



dz x dz 2 = 0.6837762984... 



and the elaborate details appear later. 

It is surprising that working with (B-A)-(D-A), (C-A)-(D-A) simultaneously 
should be so difficult. The two expressions are familiar: they are sample covariance 
coefficients 713, 723 respectively between samples of size = 3. A marginal density for 
either is found in [8], but not much else is known about off-diagonal elements of a 
Wishart matrix [9]. Of course, 713 > and 723 > if and only if the corresponding 
sample correlation coefficients P13 > and P23 > 0. A formula for a trivariate density 
for (P12, P13, P23) is outlined in jTQl [11] - evidently a sample size > 4 is presumed - 
and details still need to come together. 

4. Pinned Simplices 
A slight variation on defining a random tetrahedron in 3-space is to keep one vertex 
fixed at the origin and to select the other three vertices independently from N(0, 1) 
as before. We say that the tetrahedron is pinned. 

For the moment, let us consider pinned random Gaussian triangles ABC in 2- 
space with C = (0, 0). Note that (B - A) ■ (-A) is equal to |BJ 



1 + yg 

2 

1 + V2 



2 + v/2 y/2-v/2 5 



xl- 



2 
-1 + V2 



xl, 



l + VZ 



2-V2 A+ ^j2 + V2 B 




2 + v^ 



It follows that C falls between A and B on Q with probability 2p — 1 = l/v2. 
Equivalently, P(« > tt/2) = 1—p because cos(a) < if and only if (B — A)- (—A) < 0. 
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By symmetry, P(/3 > 7r/2) = 1 — p as well and P(7 > 7r/2) = 1/2. We deduce that 

P (pinned triangle ABC is acute) 

= l-(P(a>7r/2) + P(/3>7r/2) + P( 7 >7r/2)) 
= 1 - (2 -2p+ 1/2) = -1/2 + 1/^2 



as was to be shown. 

Let us now return to pinned random Gaussian tetrahedra in 3-space with D = 
(0, 0, 0). For brevity, we focus only on the probability that both (B — A) • (—A) > 
and (C — A) ■ {—A) > 0. Consider a simpler scenario in which a = A,b = B,c = C 
are scalars. Using Miller's [5] formulas, 



Cov 





hence 



i-i 




n=( \ I), « = -(!), " = 3, 



det(S" 1 ) = 1, -v'Z = zt + z 2 , Z'VLZ = z\ + z\, 

f(zi,z 2 ) 



x exp izi + Z2- V^yzj 



,2 i „2 
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The probability we want is given by the sextuple integral, which has value 0.834. 
but can be computed more accurately via the Fourier transform 

F(w 1 ,w 2 ) 



Mapping F(wi,w 2 ) 3 from frequency back to signal domain, we calculate 

—^- exp (zx + z 2 - V3Jz 2 + zl) dz x dz 2 = 0.8343764256... 
2V37T V / 



o o 



After a discussion of some related problems in geometric probability [12] , more 
information on F(wi,w 2 ) (also called the characteristic function corresponding to 
f(zi,z 2 )) will be given. 
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5. Random Solid Angles 
We restrict attention to pinned Gaussian random tetrahedra ABCD with D = 
(0,0,0). The dihedral angle a is the angle between normal vectors A x B, A x C 
to the triangular faces ADB, ADC respectively: 

f (A x B) ■ (A x C)\ 

a = arccos ■— — — — -- . 

\\\Ax B\\ \\Ax C\\) 

Angles (3 and 7 are defined likewise. For example, a regular tetrahedron has dihedral 
angles each equal to 

arccos(l/3) = 1.2309594173... = vr - 1.9106332362... « 70.53° « 180° - 109.47°. 

The joint density for a, (3, 7 is 



fx + y + z\ ( —x + y + z\ fx — y + z\ fx + y — 
1 cos cos cos cos 

1 V 2 y V 2 y v 2 y v 2 

71 sin(x) 2 sin(|/) 2 sin(^) 2 

iix + y + z>7i, x + y < 71 + z, y + z < 71 + x and z + x < tx + y, 
otherwise. 

As a consequence, a is uniformly distributed on [0, tx] and angles a, ft, 7 are uncor- 
rected but pairwise dependent ((«,/?) is not uniform on [0, n] x [0,7r]). 

The solid angle (or trihedral angle) at D is the area of the spherical triangle 
on the unit sphere, center D, with vertices A/ \\A\\, Bj ||S||, Gj ||C||. It is equal to 
the spherical excess a + /3 + 7 — 7T, which is between and Itx . It is also equal to 

dHIE] 

f 2 arctan(C) if C > 0, 

} 2?r + 2 arctan(C) if C < 



\A-{BxC)\ 



where 

^ = II^H ||5|| ||C|| + (A ■ B) \\C\\ + (A ■ C) \\B\\ + (B-C) \\A\\ 

It can be shown that, if a tetrahedron is acute, then each of its four solid angles are 
less than n/2, but not conversely [16J. A proposed density for the solid angle at D 
was published in 1867 [T7] : 

(x 2 - Atxx + 3tt 2 - 6) cos(x) - 6(x - 2tt) sia(x) - 2(x 2 - Anx + 3ir 2 + 3) 

16vrcos(a;/2) 4 

for < x < 27r and remained obscure until it was cited in a recent paper [18J. Details 
of the supporting geometric proof need to be carefully examined. No analytic proof 
using the joint density for a, (3, 7 has yet been found. 
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As far as is known, no analogous results are known for general Gaussian ran- 
dom tetrahedra. In particular, the sum a of the four solid angles associated with a 
tetrahedra T possesses a fascinating property [19] : 

the orthogonal projection of T onto a uniform \ 
random plane in 3-space is a triangle ) 

and it would be good to understand o more fully. As an example, the regular tetra- 
hedra has solid angles each equal to 

3 arccos(l/3) - vr = 0.5512855984..., hence a/(2n) = 0.3509593121... 

and this is the maximum such value over all equifacial tetrahedra (all faces are con- 
gruent) J2U1 ED E21 E3]- No one has studied the distribution of a when T is itself 
allowed to be random. 

While we know the mean volume of a tetrahedron [241 [25l [26l 127] with uniform 
random vertices in the unit ball (127r/715) and with uniform random vertices in the 
unit cube (3977/21600 — 7r 2 /2160), the Gaussian random scenario remains open (there 
is doubt about claims in [2"8"]). 

In 2-space, a triangle is acute if and only if its circumcenter lies inside the triangle 
(the circumcircle contains all three vertices). In 3-space, a tetrahedron is 3- well- 
centered if its circumcenter lies inside the tetrahedron; a tetrahedron is 2-well- 
centered if the circumcenter of each face lies inside the face. An acute T can fail to 
be 3-well-centered, and a 3-well-centered T can fail to be acute. However, an acute T 
must be 2-well-centered (equivalently, all its faces must be acute) but not conversely 
[TJ [29]. Many problems involving Gaussian random tetrahedra suggest themselves. 

6. Fourier Transforms 
It is not difficult to prove all the formulas we need in the "forward" direction (starting 
with / and ending with F). This is done first for the pinned case, which is easier, and 
then for the general case. Motivating the formula for the inverse Fourier transform 
of F 3 is harder. We do this for the pinned case only. 

6.1. Pinned Case: Forward Direction. Our objective is to evaluate two inte- 
grals: 



1 °? °? exp ((1 + iu)x + (1 + iv)y — \fZ\fx 1 + y 2 ) 

F(u,v) = — / / i '-dxdy, 

2tt J J \Jx l + y z 



— oo — oo 



oo oo 



G(u,v) = — -j=- l l exp ((I + iu)x + (I + iv)y — V3yx 2 + y 2 ) dxdy. 



-oo — oo 
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Let x = r cos(6), y = r sin(0), then dxdy = r dr d9 and 



1 



2tt oo 



F(u,v) = 7^—/ / exp (r (1 + iu) cos(#) + (1 + iv ) sin(0) — Vo 
o o 



dr dd 



2tt 



1 

2tt J (1 + iu) cos(6) + (l + iv) sin(0) - v^ 



de, 



2tt oo 



G(u,v) 



2^77 

1 



rexp [r 



o o 

2tt 



2^ 



7T 



(1 + i u) cos(fl) + (l + iv) sin(0) - V3 
1 



dr d6 



(1 + iu) cos(6) + (l + iv) sin(0) - V% 



T d9 



since 



Re 



ll+iu)cos(6) + {l + iv)sm{6) - Vz] = cos(0) + sin(0) - V% 

< \/2 - \/3 < 0. 



Let z = exp(i 0), then dd = —i dzj z and 



F(u,v) = — 



dz 



2tt Jc (1 + iu)\{z + \) + H + iv)±{z - \) - VS z 



1 



n Jc (l + iu){z 2 + l) -i(l + iv)(z 2 - 1) -2y/Zz 

i r 1 



nJc(v + l + iu — i)z 2 — 2y/3z — (v — 1 — iu — i 



dz 

dz, 



G(u,v) 




V5 



n Jc 



(v + 1 + i u — i)z 2 — 2y/3z — (v — 1 — iu — i) 



where C denotes the unit circle, center 0, in the complex plane. The two poles z pos , 
z QCg of each integrand are 



2^3 ± a/12 + A(v + 1 + iu - i)(v - 1 - iu - i) \/3 ± J(u - i) 2 + (v - i) 2 + 3 



2(v + l + iu-i) 



v + 1 + iu — i 
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and z neg is always inside C, z pos is always outside. For F, z neg is a pole of order 1 and 
the associated residue is 



lim 



1 



neg ( v + i + iu-i)(z-z pos ) 2 y J( u _ i y + ( v _ i y + 



multiplying by (2m)(i/'jr) completes the proof. For G, z ncg is a pole of order 2 and 
the associated residue is 



z^nog dz\(v + l + iu — i) 2 (z — z pos ) 2 



Vs 



4 [{u-i) 2 + {v-i) 2 + 3f /2 ' 

multiplying by {2m)(— 2i/(y/3ir)) completes the proof. 

6.2. General Case: Forward Direction. Our objective is to evaluate two inte- 
grals: 



F(u,v) = - 



1 



l °? °° exp f- ((1 + Aiu)x + (1 + Aiv)y- y/3y/3x' 2 -2xy + 3y 2 



dx dy, 



y/Sx 2 — 2xy + 3y 2 

— cxj — CXJ 

OO OO 

G(u,v) = — -j=- \ \ exp I - I (1 + Aiu)x + (1 + Aiv)y — v3y3x 2 — 2xy + 3y 2 ) ) dxdy. 



Let 



x 



\/2cos(0) -sin(0) 



y 



y / 2cos(^) + sin 



2v/2 v v ' ' v JJ ' 2^ 

then 3x 2 — 2x y + 3y 2 = r 2 and the Jacobian determinant is 

/ 1 \ 2 \/2cos(0) - sin(0) r(-V2sin(e)-cos(ey 
V2^/27 \/2cos(0)+sin(0) r (-\/2sin(0) + cos(^) < 



2^ 



hence dxdy = (r/(2\/2) dr d6. We obtain 

2-7T OO 



F(u,u) 



1 1 

2^2^ 



1 r 



ex P 7 







4 2\/2 



1 + Ai u) (a/2 cos(0) - sin( 



+ (l + 4^)( v / 2cos(0)+sin(0)) -2a/2-a/3 

2-7T CO 



dr d6 





2tt 



cos(0) + 2?(w + u) cos(0) - \/2i(u - v) sin(0) - y/3 
1 



dr d6 



y/2ir J Q cos(0) + 2i(u + v) cos(0) - \/2i(u - v) sin(0) - y/3 



de. 
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2tt oo 



G(u,v) 



4^2^ 



1 r 



rexp - 



o o 



4 2\/2 



— =— / / r exp , 

v 

2tt 



Vq 



TV 



(1 + 4zu) (V2cos(9) - sin(0) 

dr d6 

cos(0) + 2i(it + u) cos(6>) - >/2i(« - u) sin(0) - V3 
1 



+ (l + 4zv) (V2cos(0) + sin(0)) -2v^- v^" 

2-7T CO 



dr dd 



cos(6) + 2i(u + v) cos(6) - >/2i(« - v) sin(0) - V% 



T d9 



since 



Re 



cos(6) + 2i(u + v) cos(6>) - V2i(u - v) sin(6>) - V3 



= cos(6) - V3 
< 1 - \/3 < 0. 



Let z = exp(i 8), then d9 = —i dz/z and 
F(u,v) 



1 



dz 



V2vr Jc \{z + I) + 2z(w + u)!^ + I) - V2i(u - v ) ±-.( z - \) - v^ z 

y/2i f 1 , 

dz 



vr 7c (z 2 + 1) + 2i(u + v)(z 2 + 1) - >/2(« - v)(z 2 - 1) - 2^3* 
v^i r 1 

vr 7c (1 - y/2u + V2w + 2iu + 2iv)z 2 - 2^z + (1 + v^w - V2v + 2iu + 2iv) 



dz, 



G(u,v) 



2t 



dz 



y/6n Jc 
Si 



(z + l) + 2i(u + v)\(z + i) - >/2i(« - v)Uz -\)-y/Z 



2 



v/oV yc 

8i 



\Z6-rc Jc 



(z 2 + 1) + 2i(u + v)(z 2 + 1) - \/2(m - v)(z 2 - 1) - 2>/3^ 

2 



rrf-2 



1 - V2u + V2v + 2iu + 2iv)z 2 - 2\/Zz + (1 + v^w - v 7 ^ + 2m + 2iv) 



sdz. 



The two poles z pos , z neg of each integrand are 



2^ ± yl2 - 4(1 - V2u + v 7 ^ + 2m + 2iv)(l + y/2u - \/2v + 2%u + 2iu) 

2(1 - V2u + V2v + 2iu + 2iv) 



V3 ± y/2y/{u - i)(3« + i) + (v - i)(3v + i) + 2uv-l 
1 - \[2u + \[2v + 2m + 2m 
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and z QCg is always inside C, z pos is always outside. For F, z neg is a pole of order 1 and 
the associated residue is 

lim = ■= 

z->zneg (i _ y/2 u + ^^ + 2iu + 2iv)(i 



1 1 



2^ ^(m - j)(3 u + i) + ( v - i)(3^ + i) + 2w-l' 

multiplying by (2iii)(\f2i / it) completes the proof. For G, z ncg is a pole of order 2 and 
the associated residue is 

d ( I 

lim — - 



*-"»« cte [ (1 - v^m + V2v + 2iu + 2iv) 2 (z - z pos ) 2 j 

y/E 1 

16 [(u - i) (3m + i) + (v-i) (3v + i) + 2w - 1] 3/2 ' 

multiplying by (27ri)(— 8i/(y/Eit)) completes the proof. Clearly there is common struc- 
ture to both cases and a more encompassing theorem should be possible. 

6.3. Pinned Case: Backward Direction. The fact that G = F 3 is fairly mirac- 
ulous but completely unmotivated. Let us briefly sketch a "backward" argument for 
the pinned case only, starting with a known inverse Fourier transform [30] : 

oo . . . 

Differentiating both sides with respect to u, we obtain 

exp(-ivy)(u-i) ( r \ |l/|(«-i) 



- / «Pt-""»("-) a, = -2expfo) j f 1 f|»|V(«-i)' + 3V -J 
_i [(» - ') 2 + (« - *) 2 + 3] 3/2 V"" V( y ^(«-i)2-h3 

that is, 



7 exp(-ivy) ir 1 (|y| v /(n-z) 2 + 3) 
/ ^at; = 2 w exp ?/ . = — - 

_4 [( U - i)2 + ( V _ i)2 + 3 ]3/2 ^/( U -i) 2 + 3 

Another known inverse Fourier transform is useful now 1301: 



\ exp(—iux) . — -au 

J . (n, - i\1 J- -X 



u - % Y + 3 
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27rexp(x) / 2 2A 1 / 4 



3 V4 I,, (^ + y 2 ) 1/4 ^(V3v^T7) 

y/2nexp(x) / 2 2 \i/4 fW^V (-y/S^x 2 + y 2 



"^^.expf-v^^ + r 



Multiplying both sides by 2|y| exp(y)/(27r) 2 , we conclude that 

1 f f exp(-tux-tvy) 1 / /- r— — ; 

t — rr / / T ^dudv = — ■=- exp \x + y — v oa/x^ + y 

( 27r ) _i_i [(w-0 2 + (^-^) 2 + 3] 3/2 2V/37T P V V 

as was to be shown. 
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